
library(readstata13)

setwd("~/Dropbox/David_Paper_2/Paper Drafts/International Organization/Final Draft/Final Replication Files/Appendix Replication")

data<-read.dta13("appendix_tables_rep.dta")

un_prob<-lm(make_claim_onset_f~as.factor(year)-1,data)


plot(1816:2002,coef(un_prob),type="b",xlab="",ylab="Probability of a Claim - Dyadic",axes=F,pch=16)


axis(1,at=c(1816,1831,1848,1871,1914,1945,1991,2002))

axis(2)

grid()



### or B



data2<-read.dta13("systemic_yearly_summary2.dta")


plot(data2$year,data2$claim_rate,type="b",xlab="",ylab="Probability of a Claim",axes=F,pch=16)


axis(1,at=c(1816,1831,1848,1871,1914,1945,1991,2002))

axis(2)

grid()
